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ABSTRACT 

These notes are for the proceedings of the NATO Advanced Study Insitutue 
on Scale Invariance, Interfaces, and Non-Equilibrium Dynamics, held at the Isaac 
Newton Institute, 20-30 June 1994. The four lectures address a number of issues 
related to dynamic fluctuations of lines in non-equilibrium circumstances. The first 
two are devoted to the critical behavior of contact lines and flux lines depinning 
from impurities. It is emphasized that anisotropies in the medium lead to different 
universality classes. The importance of nonlinearities for moving lines are discussed 
On ; in the context of flux lines and polymers in the last two lectures. A dynamic form 

^ ■ birefringence is predicted for drifting polymers. 

00 

O I. DEPINNING OF A LINE IN TWO DIMENSIONS 

^j- , 

Q\ , Depinning is a non-equilibrium critical phenomenon involving an external force 

and a pinning potential. When the force is weak the system is stationary, trapped in a 
metastable state. Beyond a threshold force the (last) metastable state disappears and 
the system starts to move. While there are many macroscopic mechanical examples, 
our interest stems from condensed matter systems such as Charge Density Waves 



(CDWs) 1 , interfaces 2 , and contact lines 3 . In CDWs, the controlling parameter is the 
external voltage. A finite CDW current appears only beyond a threshold applied 
voltage. Interfaces in porous media, domain walls in random magnets, are stationary 
unless the applied force (magnetic field) is sufficiently strong. A key feature of these 
examples is that they involve the collective depinning of many degrees of freedom 
that are elastically coupled. As such these problems belong to the realm of collective 
critical phenomena, characterized by universal scaling laws. We shall introduce these 
laws and the corresponding exponents below for the depinning of a line (interface or 
contact line). 

Consider a line in two dimensions, oriented along the x direction, and fluctuating 
along a perpendicular y direction. The configuration of the line at time t is described 
by the function r(x, t). The function r is assumed to be single valued, thus excluding 
configurations with overhangs. In many cases 2 , the evolution of the curve satisfies 
an equation of the form 

*^ = F + f(x,r) + lC[r]. (1) 

The first term is a uniform applied force which is also the external control parameter. 
Fluctuations in the force due to randomness and impurities are represented by the 
second term. With the assumption that the medium is on average translationally 
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Figure 1. Geometry of the line in two dimensions. 

invariant, the average of / can be set to zero. The final term describes the elastic 
forces between different parts of the line. Short range interactions can be described 
by a gradient expansion; for example, a line tension leads to /C[r(x)] = V 2 r or 
}C[r(q)] = — q 2 r(q) for the Fourier modes. The surface of a drop of non-wetting 
liquid terminates at a contact line on a solid substrate 3 . Deformations of the contact 
line are accompanied by distortions of the liquid/gas surface. As shown by Joanny 
and de Gennes 4 , the resulting energy and forces are non-local, described by K,[r(q)] = 
— \q\r(q). More generally we shall consider the linear operator JC[r(q)] = —\q\ a r(q), 
which interpolates between the above two cases as a changes from one to two. 



Figure 2. Critical behavior of the velocity. 



When F is small, the line is trapped in one of many metastable states in which 
dr/dt = at all points. For F larger than a threshold F c , the line is depinned from 
the last metastable state, and moves with an average velocity v. On approaching the 
threshold from above, the velocity vanishes as 
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where (3 is the velocity exponent, and A is a nonuniversal amplitude. The motion 
just above threshold is not uniform, composed of rapid jumps as large segments of 
the line depin from strong pinning centers, superposed on the slower steady advance. 
These jumps have a power law distribution in size, cutoff at a correlation length £ 
which diverges at the transition as 

£~(F-F C )-". (3) 

The jumps are reminiscent of avalanches in other slowly driven systems. In fact, the 
depinning can be approached from below F c by monotonically increasing F in small 
increments, each sufficient to cause a jump to the next metastable state. The size 
and width of avalanches becomes invariant on approaching F c . For example, 

Prob(width of avalanche > £) « (4) 

where the cutoff £_ diverges as in Eq.(3). The critical line is a self-affine fractal 
whose correlations satisfy the dynamic scaling from 

([r(x, t) - r(x', t')f) = (x- x') 2C 9 ( r^L) , (5) 

defining the roughness and dynamic exponents, ( and z respectively. (Angular brack- 
ets reflect averaging over all realizations of the random force /.) The scaling function 
g goes to a constant as its argument approaches 0; £ is the wandering exponent of 
an instantaneous line profile, and z relates the average lifetime of an avalanche to its 
size by r(f) ~ £ 2 . 

Although, the underlying issues of collective depinning for CDWs and interfaces 
have been around for some time, only recently a systematic perturbative approach to 
the problem was developed. This functional renormalization group (RG) approach 
to the dynamical equations of motion was originally developed in the context of 
CDWs by Narayan and Fisher 5 (NF), and extended to interfaces by Nattermann et 
al 6 . We shall provide a brief outline of this approach starting from Eq.(l). Before 
embarking on the details of the formalism, it is useful to point out some scaling 
relations amongst the exponents which follow from underlying symmetries and non- 
renormalization conditions. 

1. As mentioned earlier, the motion of the line close to the threshold is composed 
of jumps of segments of size £. Such jumps move the interface forward by ^ over a 
time period £ z . Thus the velocity behaves as, 

v~^~\F-F c \^-0 =► = „(z-0. (6) 

2. If the elastic couplings are linear, the response of the line to a static perturbation 
s(x) is obtained simply by considering 

r £ (x,t)=r(x,t)-IC- 1 [e(x)], (7) 

where /C -1 is the inverse kernel. Since, r e satisfies Eq.(l) subject to a force F + 
e(x) + f(x,r e ), r satisfies the same equation with a force F + f(x,r — IC~ 1 [e(x)]). 
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As long as the statistical properties of the stochastic force are not modified by the 
above change in its argument, d(r)/de = 0, and 

dr e (x)\ ! m /dr c (q)\ 1 ,~ 



19 

Since it controls the macroscopic response of the line, the kernel JC cannot change 
under RG scaling. From Eqs.(5) and (3), we can read off the scaling of r(x), and the 
force SF, which using the above non-renormalization must be related by the exponent 
relation 

(+l=v- (9) 

Note that this identity depends on the statistical invariance of noise under the trans- 
formation in Eq.(7). It is satisfied as long as the correlations (f(x,r)f(x',r')) only 
depend on r - r'. The identity does not hold if these correlations also depend on the 
slope dr/dx. 

3. A scaling argument related to the Imry-Ma estimate of the lower critical dimension 
of the random field Ising model, can be used to estimate the roughness exponent 7 . 
The elastic force on a segment of length £ scales as £ < * -t7 . If fluctuations in force are 
uncorrelated in space, they scale as £ - (£ +1 )/ 2 over the area of an avalanche. Assuming 
that these two forces must be of the same order to initiate the avalanche leads to 

This last argument is not as rigorous as the previous two. Nonetheless, all three 
exponent identities can be established within the RG framework. Thus the only 
undetermined exponent is the dynamic one, z. 

A field theoretical description of the dynamics of Eq.(l) can be developed using 
the formalism of Martin, Siggia and Rose 8 (MSR): Generalizing to a d— dimensional 
interface, an auxiliary field f(x, t) is introduced to implement the equation of motion 
as a series of 5-functions. Various dynamical response and correlation functions for 
the field r(x, t) can then be generated from the functional, 

Z = J Vr{x,t)Vr{x,t)J[r]exp(S), (H) 

where 

S = i J d d xdtr(x,t){d t r - K[r\- F - f (x, r(x, £))}. (12) 

The Jacobian J[r] is introduced to ensure that the 5-functions integrate to unity. It 
does not generate any new relevant terms and will be ignored henceforth. 

The disorder-averaged generating functional Z can be evaluated by a saddle- 
point expansion around a Mean-Field (MF) solution obtained by setting /CmfH x )] = 
vt — r(x). This amounts to replacing interaction forces with Hookean springs con- 
nected to the center of mass, which moves with a velocity v. The corresponding 
equation of motion is 

±»LL = vt - r MF (t) + f[r MF (t)] + F MF (v), (13) 
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where the relationship Fmf(v) between the external force F and average velocity 
v is determined from the consistency condition (r^f(f)) = vt. The MF solution 
depends on the type of irregularity 5 : For smoothly varying random potentials, Pmf = 
3/2, whereas for cusped random potentials, (3mf = 1- Following the treatment of 
NF 5 ' 9 , we use the mean field solution for cusped potentials, anticipating jumps with 
velocity of O(l), in which case (3mf = 1- After rescaling and averaging over impurity 
configurations, we arrive at a generating functional whose low-frequency form is 

x, t)VR(x,t)exp(S), 



Z = J VR{ 

S = -J d d xdt [F — Fmf{v)] -R(x, t) 



/d d c{ duj - 

+ ^ J d d x dt dt' i?(x, t)R(x, t')C [vt - vt' + R(x, t) - R(x, t')] . 

In the above expressions, R and R are coarse-grained forms of r — vt and if, respec- 
tively. F is adjusted to satisfy the condition (R) = 0. The function C(vt) is initially 
the connected mean-field correlation function {(rMF(t)fMF(t + T ))c- 

Ignoring the /^-dependent terms in the argument of C, the action becomes Gaus- 
sian, and is invariant under a scale transformation x — > bx, t — > b a t, R — > W~ d l 2 R 1 
R -> b- a - d / 2 R, F -> 6" d / 2 F, and d — > 6 d / 2 t>. Other terms in the action, of higher 
order in R and i?, that result from the expansion of C [and other terms not explicitly 
shown in Eq.(14)], decay away at large length and time scales if d > d c = 2a. For 
d > d c , the interface is smooth (Co < 0) at long length scales, and the depinning 
exponents take the Gaussian values zq = a, vq = 2/d, (3q = 1. 

At d = d c , the action S has an infinite number of marginal terms that can be 
rearranged as a Taylor series of the marginal function C [vt — vt' + R(x, t) — -R(x, t')], 
when v — > 0. The RG is carried out by integrating over a momentum shell A/6 < 
|q| < A (we set the cutoff wave vector to A = 1 for simplicity) and all frequencies, 
followed by a scale transformation x — > bx, t — > b z t, R — > 6^i?, and — > b ~ d R, 
where b = e^. The resulting recursion relation for the linear part in the effective 
action (to all orders in perturbation theory) is 

d (F ~f MF) =( z + e)(F- F MF ) + constant, (15) 

which immediately implies (with a suitable definition of F c ) 

d(F - F c ) 



Vf{F-F c ), (16) 



di 

with the exponent identity 

y F = z + 9 = l/u . (17) 

The functional renormalization of C(u) in d = 2a — e interface dimensions, computed 
to one-loop order, gives the recursion relation, 

dC(u) r ^ n/ , > dCiu) 



(27r) d du \ du 
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where Sd is the surface area of a unit sphere in d dimensions. NF showed that all 
higher order diagrams contribute to the renormalization of C as total derivatives 
with respect to u, thus, integrating Eq.(18) at the fixed-point solution dC* jdl = 

0. together with Eqs.(9) and (17), gives ( = e/3 to all orders in e, provided that 
/ C* 7^ 0. This gives Eq.(10) for a one-dimensional interface, as argued earlier. This 
is a consequence of the fact that C{u) remains short-ranged upon renormalization, 
implying the absence of anomalous contributions to £. 

The dynamical exponent z is calculated through the renormalization of p, the 
term proportional to RdfR, which yields 

z = a -2e/9 + 0(e 2 ), (19) 

and using the exponent identity (6), 

(3 = 1 - 2e/9a + 0(e 2 ). (20) 

Nattermann et. al. 6 obtain the same results to 0(e) by directly averaging the MSR 
generating function in Eq.(ll), and expanding perturbatively around a rigidly moving 
interface. 

Numerical integration of Eq.(l) for an elastic interface 10 (a = 2) has yielded 
critical exponents ( = 0.97±0.05 and v = 1.05±0.1, in agreement with the theoretical 
result Q = v = 1. The velocity exponent (3 = 0.24 ±0.1 is also consistent with the 
one-loop theoretical result 1/3; however, a logarithmic dependence v ~ 1/ ln(F — F c ), 
which corresponds to (3 = 0, also describes the numerical data well. In contrast, 
experiments and various discrete models of interface growth have resulted in scaling 
behaviors that differ from system to system. A number of different experiments 
on fluid invasion in porous media 11 give roughness exponents of around 0.8, while 
imbibition experiments 12 ' 13 have resulted in ( « 0.6. A discrete model studied by 
Leschhorn 14 , motivated by Eq.(l) with a = 2, gives a roughness exponent of 1.25 
at threshold. Since the expansion leading to Eq.(l) breaks down when £ approaches 
one, it is not clear how to reconcile the results of Leschhhorn's numerical work 14 with 
the coarse-grained description of the RG calculation, especially since any model with 
( > 1 cannot have a coarse grained description based on gradient expansions. 

Amaral, Barabasi, and Stanley (ABS) 15 recently pointed out that various models 
of interface depinning in 1+1 dimensions fall into two distinct classes, depending on 
the tilt dependence of the interface velocity: 

1. For models like the random field Ising Model 16 , and some Solid On Solid models, 
the computed exponents are consistent with the exponents given by the RG analysis. 
It has been suggested 14 , however, that the roughness exponent is systematically larger 
than e/3, casting doubt on the exactness of the RG result. 

2. A number of different models, based on directed percolation (DP) 17 ' 12 give a 
different roughness exponent, ( ~ 0.63. In these models, pinning sites are randomly 
distributed with a probability p, which is linearly related to the force F. The interface 
is stopped by the boundary of a DP cluster of pinning sites. The critical exponents at 
depinning can then be related to the longitudinal and transverse correlation length 
exponents v\\ « 1.70 and v±_ w 1.07 of DP. In particular, Q = v\\/v±_ ~ 0.63, and 
j3 = zv|| — v± ~ 0.63, in agreement with experiments. 
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The main difference of these models can be understood in terms of the depen- 
dence of the threshold force F c to the orientation. To include the possible dependence 
of the line mobility on its slope, d x r, we can generalize the equation of motion to 

d t r = Kd 2 x r + Kd x r + -(d x r) 2 + F + f(x, r). (21) 

The isotropic depinning studied by RG corresponds to k = A = 0. In models of 
depinning by directed percolation studied so far 17 ' 12 there is a dependence of F c on 
slope, making a nonzero A unavoidable. The nonlinearity is relevant, accounting for 
the different universality class. Eq.(21) with k = 0, motivated in a different fashion, 
has been studied by Stepanow 18 . The exponents obtained approximately by a one 
loop expansion, £ ps 0.8615, z = 1, and yp ~ 0.852 are reasonably close to those 
of directed percolation. The presence of anisotropy in depinning actually suggests a 
third possibility: 

3. When the line is depinning along a (tilted) direction of lower symmetry, even 
more relevant terms like nd x r will be present in the equation of motion. This new 
universality class is possibly controlled by "tilted" DP clusters 19 , for which £ = 1/2. 

For the case of the contact line (CL) (a = 1), these anisotropies are irrelevant, 
but there are other concerns related to the details of the driving force: In most 
experiments, the velocity of the CL is controlled rather than the external force. The 
effect of this can be numerically investigated by replacing the external force F in 
Eq.(l) with 

F' = v- J ^f(x>,r(x',t)), (22) 

and looking at the time average of F' as a function of v. (F' is chosen such that 
J dxd t r(x) = vL.) Even though the critical behavior for both ways of driving may 
be the same for an infinitely large system, there is a system size dependent region 
near the depinning threshold where the behavior changes drastically. Preliminary 
findings on an elastic line suggest that in this region, the velocity exponent (3 be- 
comes considerably larger than one, in marked contrast with the constant force case. 
This can be qualitatively understood as follows: For a system of finite size, when a 
constant driving force is applied, the average velocity drops to zero as soon as tempo- 
ral fluctuations of the instantaneous velocity are comparable with the time-averaged 
velocity. This is because the time average is then completely dominated by configu- 
rations for which the interface is pinned. Thus, the pinning transition becomes truly 
second order only in the large system limit: The velocity jumps to zero from a finite 
value in a finite system. In contrast, for constant velocity driving, no configuration 
has more weight than any other, since the interface is constrained to move past any 
obstacles by suitably increasing the applied external force, and decreasing it when 
passing through weakly pinning regions. Thus, in the region where a force-driven 
interface is pinned, the velocity-driven interface will experience fluctuations in the 
external force comparable to the average force itself. This average force as a function 
of velocity has an effective velocity exponent much larger than one. This distinction 
may partially explain the large velocity exponent found in a recent CL experiment 20 , 
where the interface was velocity-driven. In addition to this, gravity imposes a finite 
wavelength cutoff on the roughening of the CL, which may complicate the analysis 
of experimental results. 
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II. DEPINNING OF A LINE IN THREE DIMENSIONS 

The pinning of flux lines (FLs) in Type-II superconductors is of fundamental 
importance to many technological applications that require large critical currents 21 . 
Upon application of an external current density J, the FL becomes subject to a 
Lorentz force per unit length 

F = ^Jxt, (23) 
c 

where <po is the flux quantum, and t is the unit tangent vector along the FL, which 
points along the local magnetic field. The motion of FLs due to the Lorentz force 
causes undesirable dissipation of supercurrents. Major increases in the critical current 
density J c of a sample are achieved when the FLs are pinned to impurities. 

Recent numerical simulations have concentrated on the low temperature behav- 
ior of a single FL near depinning 22 ' 10 ' 23 , mostly ignoring fluctuations transverse to 
the plane defined by the magnetic field and the Lorentz force. Common signatures 
of the depinning transition from J < J c to J > J c include a broad band (f~ a type) 
voltage noise spectrum, and self-similar fluctuations of the FL profile. 




Figure 3. Geometry of the line in three dimensions. 



The configuration of the FL at time t is now described by the vector function 
r(x,t), where x is along the magnetic field B, and the unit vector ey is along the 
Lorentz force F. (See Fig. 3) The major difference of the FL from the line in two 
dimensions is that the position, r(x, t), is now a 2-dimensional vector instead of a 
scalar; fluctuating along both ey and e_i_ directions. Point impurities are modeled by a 
random potential V(x, r), with zero mean and short-range correlations. The simplest 
possible Langevin equation for the FL, consistent with local, dissipative dynamics, is 

p-^ = c» 2 r + f(x,r(x,t))+F, (24) 
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where p is the inverse mobility of the FL, and f = — V r V . The potential V(x, r) need 
not be isotropic. For example, in a single crystal of ceramic superconductors with the 
field along the oxide planes, it will be easier to move the FL along the planes. This 
leads to a pinning threshold that depends on the orientation of the force. Anisotropy 
also modifies the line tension, and the elastic term in Eq.(24) is in general multiplied 
by a non-diagonal matrix K a p. The random force f (x, r), can be taken to have zero 
mean with correlations 

(f a (x, r)/ 7 (x', r')> = S(x - x')A ai (r - r'). (25) 

We shall focus mostly on the isotropic case, with A Q7 (r — r') = <5 Q7 A(|r — r'|), where 
A is a function that decays rapidly for large values of its argument. 

In addition to the exponents defined in the first lecture, (/?, v, £ — > z ~~ z \\), 
there are two additional critical exponents that describe fluctuations transverse to 
the overall motion of the FL slightly above depinning. At length scales up to £, the 
correlated fluctuations satisfy the dynamic scaling form, 

([r ll (x,t)-r ll (x',t')] 2 ) = \x-x'\ 2 ^g ll 
([rx(x, t) - r ± (x\ t')f) = \x- x'\ 2 ^g ± 

where (j_ and z± are the transverse roughness and dynamic exponents. One conse- 
quence of transverse fluctuations is that a "no passing" rule 24 , applicable to CDWs 
and interfaces, does not apply to FLs. It is possible to have coexistence of moving 
and stationary FLs in particular realizations of the random potential. The effects 
of transverse fluctuations r± for large driving forces, when the impurities act as 
white noise, will be discussed later. At this point, we would like to know how these 
transverse fluctuations scale near the depinning transition, and whether or not they 
influence the critical dynamics of longitudinal fluctuations near threshold. 

The answer to the second question is obtained by the following qualitative argu- 
ment: Consider Eq.(24) for a particular realization of randomness f(x,r). Assuming 
that portions of the FL always move in the forward direction 25 , there is a unique 
point r±(x,r\\) that is visited by the line for given coordinates (x,r\\). We con- 
struct a new force field /' on a two dimensional space (x,r\\) through f'(x,r\\) = 
/|| (x,r\\, r±(x, r||)). It is then clear that the dynamics of the longitudinal component 
7"|| (x, t) in a given force field f (x, r) is identical to the dynamics of ry (x, t) in a force 
field f'(x,r\\), with r± set to zero. It is quite plausible that, after averaging over 
all f, the correlations in /' will also be short-ranged, albeit different from those of 
f . Thus, the scaling of longitudinal fluctuations of the depinning FL will not change 
upon taking into account transverse fluctuations. However, the question of how these 
transverse fluctuations scale still remains. 

Certain statistical symmetries of the system restrict the form of response and cor- 
relation functions. For example, Eq.(24) has statistical space- and time-translational 
invariance, which enables us to work in Fourier space, i.e. (x,t) — > (q,u). For an 
isotropic medium, F and v are parallel to each other, i.e., v(F) = v(F)F, where F 
is the unit vector along F. Furthermore, all expectation values involving odd powers 
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of a transverse component are identically zero due to the statistical invariance un- 
der the transformation r± — > — r±. Thus, linear response and two-point correlation 
functions are diagonal. The introduced critical exponents are then related through 
scaling identities. These can be derived from the linear response to an infinitesimal 
external force field e(q,u>), 

/ dr a (q,u) \ _ 

in the (q, uj) — > (0, 0) limit. Eq.(24) is statistically invariant under the transformation 
F — > F + e{q), v(q,u) — > v(q,u) + q~ 2 e(q). Thus, the static linear response has the 
form x\\(Qi 00 = 0) = x±( a i UJ = 0) = q~ 2 . Since s\\ scales like the applied force, the 
form of the linear response at the correlation length £ gives an exponent identity 
similar to Eq.(9): 

C|| + = 2. (28) 

Considering the transverse linear response seems to imply (± = Q. However, the 
static part of the transverse linear response is irrelevant at the critical RG fixed 
point, since z± > z\\ , as shown below. When a slowly varying uniform external force 
e(t) is applied, the FL responds as if the instantaneous external force F + e is a 
constant, acquiring an average velocity, 

f)i) 

(d t r a ) = v a (F + e)^ v a (F) + ^£ 7 . (29) 

Substituting dv\\/dF\\ = dv/dF and dv±/dF± = v/F, and Fourier transforming, 
gives 

1 

X\\(q = 0,u)) = — ■ /, /jtt\-i i 7v 2^ ' 

-iu(dv/dF) 1 + 0{u 2 ) ^ 



X±(q = 0,u>) = 



-iiviv/F)- 1 + 0(lv 2 )' 



Combining these with the static response, we see that the characteristic relaxation 
times of fluctuations with wavelength £ are 

which, using Eq.(28) , yield the scaling relations 

/? = (*H-C||K (32) 

z± = z\\ + 1/v. 

We already see that the dynamic relaxation of transverse fluctuations is much slower 
than longitudinal ones. All critical exponents can be calculated from and z\\, 

by using Eqs.(28) and (32). 
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Equation (24) can again be analyzed using the MSR formalism. The long wave- 
length, low frequency behavior, for isotropic random potentials, is described by the 
effective action 



S = - J dtd d x[F - F M f(v)]R\\(x, t) 

d d q doj - . . . . . . 2 \ 



(27r) d 2tt 



i2ii(-q,-o;)fl||(q,u;)(-iu;p + g ) 



d d q dw~ . . _ . . / F c 2 
-Rj_(-q, -u) ■ Rj_(q, <jJ) -iu) h q 



(27r) d 2/T \ v 

+ / d d xdtdt' R 1 {^t)R 1 {^t')C 1 {v(t-t') + J R||(x,£) -i?||(x,t')) . 

(33) 

All terms in 5 involving longitudinal fluctuations are identical to the two-dimensional 
case, thus we obtain the same critical exponents for longitudinal fluctuations, i.e., 
C|| = e/3, Z|| = 2 — 2e/9 + 0(e 2 ). The renormalization of transverse temporal force- 
force correlations C±(u) yields an additional recursion relation 

^M = [£ + 2 ^ + 2 (, || -2)]C ± ( ll ) +C|| «^M 

(34) 

In addition, the form of the transverse dynamic linear response given in Eq.(30) 
implies the nonrenormalization of the term proportional to R±dtR±(F c /v), which, 
along with the renormalization of Cj_(0), gives a transverse roughness exponent (± = 
C\\ — d/2, to all orders in perturbation theory. For the FL (e = 3), the critical 
exponents are then given by 

C|| = 1, z,, « 4/3, v = 1, 
/? ~ 1/3, C± = 1/2, zl«7/3. 

Numerical integrations of Eq.(24) 26 that test the scaling forms and exponents 
predicted by Eqs.(2) and (26) are in agreement with RG results: A fit for the velocity 
exponent gives (3 = 0.3 ± 0.1, although a logarithmic fit {j3 = 0) cannot be ruled out, 
as seen in Fig. 4. The roughness exponents (see Fig. 5) fit the scaling form well, with 10 
C|| = 0.94 ± 0.05, and ( ± = 0.50 ± 0.02. 

The potential pinning the FL in a single superconducting crystal is likely to 
be highly anisotropic. For example, consider a magnetic field parallel to the copper 
oxide planes of a ceramic superconductor. The threshold force then depends on its 
orientation, with depinning easiest along the copper oxide planes. In general, the 
average velocity may depend on the orientations of the external force and the FL. 
The most general gradient expansion for the equation of motion is then, 

dr 1 

= fiapFp + K a0 d x ri3 + K a pdlrp + -K,p~ i d x r p d x r 1 + f a (x, r(x, £)) H , (36) 
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• Numerical data 
— Fit to (F-F c ) f , p s 0.31 
-- Fit to |ln(F/F c -1)| a , a = 0.64 



0.0 



1.00 1.10 
External Force 



Figure 4. A plot of average velocity versus external force for a sys- 
tem of 2048 points. Statistical errors are smaller than symbol sizes. 
Both fits have three adjustable parameters: The threshold force, the 
exponent, and an overall multiplicative constant. 




Figure 5. A plot of equal time correlation functions versus sepa- 
ration, for the system shown in Fig. 4, at F = 0.95. The observed 
roughness exponents very closely follow the theoretical predictions 
of C|| =1, Cl = 0.5, which are shown as solid lines for comparison. 



with 

(f a (x, r)f p (x', r')) = S(x - x')C a p{v - r'). (37) 

Depending on the presence or absence of various terms allowed by the symmetries of 
the system, the above set of equations encompasses many distinct universality classes. 
For example, consider the situation where v depends on F, but not on the orientation 
of the line. Eqs.(27) and (30) have to be modified, since v and F are no longer 
parallel (except along the axes with r — > — r symmetry), and the linear response 
function is not diagonal. The RG analysis is more cumbersome: For depinning along 
a non-symmetric direction, the longitudinal exponents are not modified (in agreement 
with the argument presented earlier), while the transverse fluctuations are further 
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suppressed to (j_ = — 2 (equal to zero for Q = l) 27 . Relaxation of transverse 
modes are still characterized by z± = z» + and the exponent identity (28) also 
holds. Surprisingly, the exponents for depinning along axes of reflection symmetry 
are the same as the isotropic case. If the velocity also depends on the tilt, there 
will be additional relevant terms in the MSR partition function, which invalidate the 
arguments leading to Eqs.(28)-(32). The analogy to FLs in a planes suggests that the 
longitudinal exponents for d = 1 are controlled by DP clusters ' 12 , with £y m 0.63. 
Since no perturbative fixed point is present in this case, it is not clear how to explore 
the behavior of transverse fluctuations systematically. 

III. NONLINEAR DYNAMICS OF MOVING LINES 

We have so far investigated the dynamics of a line near the depinning transition. 
Now, we would like to consider its behavior in a different regime, when the external 
driving force is large, and the impurities appear as weak barriers that deflect portions 
of the line without impeding its overall drift. In such non-equilibrium systems, one 
can regard the evolution equations as more fundamental, and proceed by constructing 
the most general equations consistent with the symmetries and conservation laws of 
the situation under study 28 . Even in a system with isotropic randomness, which we 
will discuss here, the average drift velocity, v, breaks the symmetry between forward 
and backward motions, and allows introduction of nonlinearities in the equations of 
motion 29 ' 28 . 

Let us first concentrate on an interface in two dimensions. (Fig.l.) By con- 
tracting up to two spatial derivatives of r, and keeping terms that are relevant, one 
obtains the Kardar-Parisi-Zhang 30 (KPZ) equation, 

d t r(x, t) = fiF + Kd 2 x r(x, t) + ^ [d x r(x, t)f + f(x, t), (38) 

with random force correlations 

(f(x, t)f(x', f )) = 2T5(x - x')5(t - t'). (39) 

For a moving line, the term proportional to the external force can be absorbed without 
loss of generality by considering a suitable Galilean transformation, r — > r — at, to a 
moving frame. A large number of stochastic nonequilibrium growth models, like the 
Eden Model and various ballistic deposition models are known to be well described, 
at large length scales and times, by this equation, which is intimately related to 
several other problems. For example, the transformation v(x,t) = —\d x r(x,t) maps 
Eq.(38) to the randomly stirred Burgers' equation for fluid flow 31 ' 32 , 

d t v + vd x v = Kd 2 x v-\d x f(x,t). (40) 

The correlations of the line profile still satisfy the dynamic scaling form in Eq.(5), 
nevertheless with different scaling exponents £, z and scaling function g. This self- 
affine scaling is not critical, i.e., not obtained by fine tuning an external parameter 
like the force, and is quite different in nature than the critical scaling of the line near 
the depinning transition, which ceases beyond the correlation length scale £. 

Two important nonperturbative properties of Eq.(38) help us determine these 
exponents exactly in 1+1 dimensions: 
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1. Galilean Invariance (GI): Eq.(38) is statistically invariant under the infinitesimal 
reparametrization 

r' = r + ex , x' = x + Xet , t' = t, (41) 

provided that the random force / does not have temporal correlations 33 . Since the 
parameter A appears both in the transformation and Eq.(38), it is not renormalized 
under any RG procedure that preserves this invariance. This implies the exponent 
identity 32 ' 33 

C + z = 2. (42) 

2. Fluctuation-Dissipation (FD) Theorem: Eqs.(38) and (39) lead to a Fokker- 
Planck equation for the evolution of the joint probability V [r(x)], 

f ( SV 5 2 V \ 

d ' v =j dx {m dtr+T w)- (43) 

It is easy to check that V has a stationary solution 

P = 6XP ("lf / dx ^ 2 )- ( 44 ) 

If V converges to this solution, the long-time behavior of the correlation functions 
in Eq.(5) can be directly read off Eq.(44), giving Q = 1/2. 

Combining these two results, the roughness and dynamic exponents are exactly 
determined for the line in two dimensions as 

C = 1/2 , z = 3/2. (45) 

Many direct numerical simulations and discrete growth models have verified these 
exponents to a very good accuracy. Exact exponents are not known for interfaces 
in higher dimensions, since the FD property is only valid in two dimensions. These 
results have been summarized in a number of recent reviews 34 ' 35 . 

Let us now turn to the case of a line in three dimensions (Fig. 3). Fluctuations 
of the line can be indicated by a a two dimensional vector r. Even in an isotropic 
medium, the drift velocity v breaks the isotropy in r by selecting a direction. A 
gradient expansion up to second order for the equation of motion gives 36 



d x r l3 d x r 1 



(46) 



d t r a = [if i<5 a/3 + K 2 v a vp] dlrp 

+ [Ai(5 Q/ 3t> 7 + 8 al Vf3) + \2V ce 5 f 3 1 + \3VaV(3V^} — 1 + f a 

with random force correlations 

(f a (x, t)U(x', t')) = 2{T 1 5 otl3 + T 2 v a v p ]5(x - x')5{t - t'). (47) 

Higher order nonlinear ities can be similarly constructed but are in fact irrelevant. In 
terms of components parallel and perpendicular to the velocity, the equations are 

dtn =K\\dl n + ^{d xn ) 2 + ^{d x r^) 2 + f^t) 

1 1 , (48) 

d t r± = K±dlr± + \±d x r\\d x r± + f±(x,t) 
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with 

(f ll (x,t)f ll (x',t , ))=2T l] 6(x-x')8(t-t') 

(f±(x,t)f ± (x',t')) =2T ± S(x-x')S(t-t') 
The noise-averaged correlations have a dynamic scaling form like Eq.(26), 

(ir ll (x,t)-r ll (x',t')] 2 ) = \x-x'\ 2 ^g ll ( 
([r ± (x,t)-r ± (x',t')f) = \x-xf^g ± ( |j!_"J|L ) • 





t-t' 




\x 


— x'\ 


z \\ 



(49) 



(50) 




Figure 6. A projection of RG flows in the parameter space, for n = 1 
transverse components. 



In the absence of nonlinearities (Ay = A x = Aj_ = 0), Eqs.(48) can easily be 
solved to give £|| = (j_ = 1/2 and z\\ = z± = 2. Simple dimensional counting indicates 
that all three nonlinear terms are relevant and may modify the exponents in Eq.(50). 
Studies of related stochastic equations 37 ' 38 indicate that interesting dynamic phase 
diagrams may emerge from the competition between nonlinearities. Let us assume 
that Ay is positive and finite (its sign can be changed by ry — > —r\\), and focus on 
the dependence of the scaling exponents on the ratios X±/X\\ and A x /A||, as depicted 
in Fig. 6. (It is more convenient to set the vertical axis to \ X K\\T_\_/\\\K±T\\.) 

The properties discussed for the KPZ equation can be extended to this higher 
dimensional case: 

1. Galilean Invariance (GI): Consider the infinitesimal reparametrization 

x' = x + Anet , t' = t , 

, " , ( 51 ) 

r\\ = r\\ + ex , r± = r± . 
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Eqs.(48) are invariant under this transformation provided that Ay = Aj_. Thus along 
this line in Fig. 6 there is GI, which once more implies the exponent identity 



C||+*||=2. 



(52) 



2. Fluctuation-Dissipation (FD) Condition: The Fokker-Planck equation for the 
evolution of the joint probability V [r\\ (x), r±(x)~\ has a stationary solution 



Vq oc exp 



dx 



|i (W + £ (W ]) 



(53) 



provided that X X K\\T± = \±K±T\\. Thus for this special choice of parameters, 
depicted by a starred line in Fig. 6, if V converges to this solution, the long-time 
behavior of the correlation functions in Eq.(50) can be directly read off Eq.(53), 
giving C|| = C± = V 2 - 

3. The Cole-Hopf ( CH) Transformation is an important method for the exact study 
of solutions of the one component nonlinear diffusion equation 31 . Here we generalize 
this transformation to the complex plane by defining, for A x < 0, 



ty(x, t) = exp 



Ayry (x, t) + iy/-\\\\ x r±(x, t) 



2K 



(54) 



The linear diffusion equation 

d t ^ = KdZ* + n(x,t)V, 

then leads to Eqs.(48) if K\\ = K±_ = K and Ay = Aj_. [Here Re(p) = X\\f\\/2K 
and Im( / u) = ^—X\\X x fj_/2K.] This transformation enables an exact solution of the 
deterministic equation, and further allows us to write the solution to the stochastic 
equation in the form of a path integral 



(x,t) r rt 

Vx(t) exp < — / dr 
(o,o) I Jo 



X 

2K 



+ fj,(x,r) 



(55) 



Eq.(55) has been extensively studied in connection with quantum tunneling in a dis- 
ordered medium 39 , with \1/ representing the wave function. In particular, results for 
the tunneling probability |\1/| 2 suggest z\\ =3/2 and £|| = 1/2. The transverse fluctu- 
ations correspond to the phase in the quantum problem which is not an observable. 
Hence this mapping does not provide any information on (± and z± which are in fact 
observable for the moving line. 

At the point Aj_ = A x = 0, r\\ and r±_ decouple, and z±_ = 2 while z\\ = 3/2. 
However, in general z\\ = z± = z unless the effective Xj_ is zero. For example at 
the intersection of the subspaces with GI and FD the exponents z\\ = z± = 3/2 
are obtained from the exponent identities. Dynamic RG recursion relations can be 
computed to one-loop order 36,40 , by standard methods of momentum-shell dynamic 
RG 32 ' 33 . 
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The renormalization of the seven parameters in Eqs.(48), generalized to n trans- 
verse directions, give the recursion relations 



de 

dK ± 



K± 



1 \±XxT ± 



n 1 A S T H 
9 -I - 4- n- 

71 4i^| 7T 4if||if2 

1 A ± ((A X T ± /K ± ) + (A ± T||/K||)) 



2-2 + 



7T 2K ± (K ± + K ll ) 

1 #x Ax((A x Tx/#x) - (A ± T||/K||)) 



7T Kj_ + #11 



#x(#X+#||) 



dZj| 

cffx 
~~df 

de 

d\ ± 

If 

dXx 

de 



z - 2C|| - 1 + 



20 



1 + 



1^1 

7T 4# 3 
1 



+ n 



1 A x Tf 



7T AK\ 



7T + 



A ll [C||+«-2], 



= A ± 



A, 



C|| + ^ - 2 



1 A,, - Ax 



7T (# ± 
2C± - C|| + * - 2 + 



((A x Tx/#x)-(AxT||/#||)) 



+ K \\) 2 
1 A||#x-Ax#|| 

7T#x#||(#X +#||) 



((A x Tx/#x)-(AxT||/#||)) 



(56) 

The projections of the RG flows on the two parameter subspace shown in Fig. 6 
are indicated by trajectories. They naturally satisfy the constraints imposed by 
the non-perturbative results: the subspace of GI is closed under RG, while the FD 
condition appears as a fixed line. The RG flows, and the corresponding exponents, 
are different in each quadrant of Fig. 6, which implies that the scaling behavior is 
determined by the relative signs of the three nonlinearities. This was confirmed by 
numerical integrations 36 ' 40 of Eqs.(48), performed for different sets of parameters. A 
summary of the computed exponents are given in Table I. 

The analysis of analytical and numerical results can be summarized as follows: 

AxA x > : In this region, the scaling behavior is understood best. The RG 
flows terminate on the fixed line where FD conditions apply, hence £|| = (j_ = 1/2. 
All along this line, the one loop RG exponent is z = 3/2. These results are consistent 
with the numerical simulations. The measured exponents rapidly converge to these 
values, except when Ax or A x are small. 

A x =0: In this case the equation for r\\ is the KPZ equation (38), thus (\\ = 1/2 
and z\\ = 3/2. The fluctuations in ry act as a strong (multiplicative and correlated) 
noise on r±. The one-loop RG yields the exponents z± = 3/2, (x = 0.75 for Ax > 0, 
while a negative Ax scales to suggesting z± > z\\ . Simulations are consistent with 
the RG calculations for Ax > 0, yielding (x = 0.72, surprisingly close to the one-loop 
RG value. For Ax < 0, simulations indicate z± ~ 2 and (± ~ 2/3 along with the 
expected values for the longitudinal exponents. 

Ax = 0: The transverse fluctuations satisfy a simple diffusion equation with 
Cx = 1/2 and z± = 2. Through the term A x (<9 x rx) 2 /2, these fluctuations act as 
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a correlated noise 33 for the longitudinal mode. A naive application of the results 
of this reference 33 give £|| = 2/3 and z\\ = 4/3. Quite surprisingly, simulations 
indicate different behavior depending on the sign of A x . For A x < 0, z\\ ~ 3/2 and 
Cy ~ 1/2 whereas for A x > 0, longitudinal fluctuations are much stronger, resulting in 
z\\ ^ 1.18 and C|| ~ 0.84. Actually, £|| increases steadily with system size, suggesting a 
breakdown of dynamic scaling, due to a change of sign in Aj_A x • This dependence on 
the sign of A x may reflect the fundamental difference between behavior in quadrants 
II and IV of Fig.6. 

TABLE I. Numerical estimates of the scaling exponents, for various val- 
ues of model parameters for n = 1. In all cases, K\\ = K± = 1 and 
Tj| = Tj_ = 0.01, unless indicated otherwise. Typical error bars are ±0.05 
for (, ±0.1 for z/(. Entries in brackets are theoretical results. Exact values 
are given in fractional form. 
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20 
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0.50 4.0 
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< z±) 


(1/2) (4) 


20 


-20 





0.55 


2.9 


0.51 4.0 








(-11 


< z ± ) 


(1/2) (4) 



Aj^ < and A x > 0: The analysis of this region (II) is the most difficult in that 
the RG flows do not converge upon a finite fixed point and Aj_ — > 0, which may 
signal the breakdown of dynamic scaling. Simulations indicate strong longitudinal 
fluctuations that lead to instabilities in the discrete integration scheme, excluding 
the possibility of measuring the exponents reliably. 

Aj^ > and A x < 0: The projected RG flows in this quadrant (IV) converge to 
the point Aj_/A|| = 1 and A x T^if||/A||Tj|.ftT_|_ = —1. This is actually not a fixed point, 
as K\\ and K± scale to infinity. The applicability of the CH transformation to this 
point implies z\\ = 3/2 and (\\ = 1/2. Since A^ is finite, z± = z\\ = 3/2 is expected, 
but this does not give any information on Simulations indicate strong transverse 
fluctuations and suffer from difficulties similar to those in region II. 

Eqs.(48) are the simplest nonlinear, local, and dissipative equations that govern 
the fluctuations of a moving line in a random medium. They can be easily generalized 
to describe the time evolution of a manifold with arbitrary internal (x e R d ) and 
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external (r e R n+1 ) dimensions, and to the motion of curves that are not necessarily 
stretched in a particular direction. Since the derivation only involves general sym- 
metry arguments, the given results are widely applicable to a number of seemingly 
unrelated systems. We will discuss one application to drifting polymers in more detail 
in the next lecture, explicitly demonstrating the origin of the nonlinear terms start- 
ing from more fundamental hydrodynamic equations. A simple model of crack front 
propagation in three dimensions 41 also arrives at Eqs.(48), implying the self-affine 
structure of the crack surface after the front has passed. 

IV. NONLINEAR RELAXATION OF DRIFTING POLYMERS 

The dynamics of polymers in fluids is of much theoretical interest and has been 
extensively studied 42 ' 43 . The combination of polymer flexibility, interactions, and 
hydrodynamics make a first principles approach to the problem quite difficult. There 
are, however, a number of phenomenological studies that describe various aspects of 
this problem 44 . 




Figure 7. The configuration of a polymer. 

One of the simplest is the Rouse model 45 : The configuration of the polymer at 
time t is described by a vector R(x,t), where x G [0, N] is a continuous variable 
replacing the discrete monomer index (see Fig. 7). 

Ignoring inertial effects, the relaxation of the polymer in a viscous medium is 
approximated by 

d t R{x, t) = fiF(R(x, t)) = Kd 2 x ~R(x, t) + rj(x, t), (57) 

where \x is the mobility. The force F has a contribution from interactions with near 
neighbors that are treated as springs. Steric and other interactions are ignored. The 
effect of the medium is represented by the random forces r\ with zero mean. The 
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Rouse model is a linear Langevin equation that is easily solved. It predicts that the 
mean square radius of gyration, K = (|R— (R) | 2 ), is proportional to the polymer size 
N, and the largest relaxation times scale as the fourth power of the wave number, 
(i.e., in dynamic light scattering experiments, the half width at half maximum of 
the scattering amplitude scales as the fourth power of the scattering wave vector 
q). These results can be summarized as R g ~ N v and T(q) ~ q z , where v and z 
are called the swelling and dynamic exponents, respectively 46 . Thus, for the Rouse 
Model, v = 1/2 and z = 4. 

The Rouse model ignores hydrodynamic interactions mediated by the fluid. 
These effects were originally considered by Kirkwood and Risemann 47 and later on 
by Zimm 48 . The basic idea is that the motion of each monomer modifies the flow field 
at large distances. Consequently, each monomer experiences an additional velocity 

s H d t R(*,t) = j_ r dX ' F{x ' )r -' tgM^ - r dx > \, din, 



8nr] s J |r xx '| 3 J \x 

(58) 

where v xx > = R(x) — R(x') and the final approximation is obtained by replacing 
the actual distance between two monomers by their average value. The modified 
equation is still linear in R and easily solved. The main result is the speeding up of 
the relaxation dynamics as the exponent z changes from 4 to 3. Most experiments 
on polymer dynamics 49 indeed measure exponents close to 3. Rouse dynamics is still 
important in other circumstances, such as diffusion of a polymer in a solid matrix, 
stress and viscoelasticity in concentrated polymer solutions, and is also applicable to 
relaxation times in Monte Carlo simulations. 

Since both of these models are linear, the dynamics remains invariant in the 
center of mass coordinates upon the application of a uniform external force. Hence 
the results for a drifting polymer are identical to a stationary one. This conclusion 
is in fact not correct due to the hydrodynamic interactions. For example, consider 
a rodlike conformation of the polymer with monomer length bo where d x R a = bot a 
everywhere on the polymer, so that the elastic (Rouse) force vanishes. If a uniform 
force E per monomer acts on this rod, the velocity of the rod can be solved using 
Kirkwood Theory, and the result is 42 

v=t^E-[I + tt]. ( 59 ) 

In the above equation, n s is the solvent viscosity, t is the unit tangent vector, k = 
2b/boN is the ratio of the width b to the half length boN/2 of the polymer. A 
more detailed calculation of the velocity in the more general case of an arbitrarily 
shaped slender body by Khayat and Cox 50 shows that nonlocal contributions to the 
hydrodynamic force, which depend on the whole shape of the polymer rather than 
the local orientation, are 0(1/(\q.k) 2 ). Therefore, corrections to Eq.(59) are small 
when N ^>b/b . 

Incorporating this tilt dependence of polymer mobility requires adding terms 
nonlinear in the tilt, <9 x r, to a local equation of motion. Since the overall force 
(or velocity) is the only vector breaking the isotropy of the fluid, the structure of 
these nonlinear terms must be identical to eq.(46). Thus in terms of the fluctuations 
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parallel and perpendicular to the average drift, we again recover the equations, 



where {_L i} refers to the 2 transverse coordinates of the monomer positions. The 
noise is assumed to be white and gaussian but need not be isotropic, i.e. 



At zero average velocity, the system becomes isotropic and the equations of motion 
must coincide with the Rouse model. Therefore, {Ay, A x , Aj_, U, K\\ — K±, Tj| — T±} 
are all proportional to E for small forces. The relevance of these nonlinear terms are 
determined by the dimensionless scaling variable 



where U* is a characteristic microscopic velocity associated with monomer motion 
and is roughly 10-20 m/s for polystyrene in benzene. The variable y is proportional 
to another dimensionless parameter, the Reynolds number Re, which determines the 
breakdown of hydrodynamic equations and onset of turbulence. However, typically 
Re <C y, and the hydrodynamic equations are valid for moderately large y. Eqs. (60) 
describe the static and dynamical scaling properties of the nonlinear and anisotropic 
regime when U > U*N~ 1 / 2 . 

Eq.(60) is just a slight variation from (48), with two transverse components 
instead of one. Thus, the results discussed in the previous lecture apply. A more 
detailed calculation of the nonlinear terms from hydrodynamics 51 shows that all three 
nonlinearities are positive for small driving forces. In this case, the asymptotic scaling 
exponents are isotropic, with v = 1/2 and z = 3. However, the fixed points of the 
RG transformation are in general anisotropic, which implies a kinetically induced 
form birefringence in the absence of external velocity gradients. This is in marked 
contrast with standard theories of polymer dynamics where a uniform driving force 
has essentially no effect on the internal modes of the polymer. 

When one of the nonlinearities approaches to zero, the swelling exponents may 
become anisotropic and the polymer elongates or compresses along the longitudinal 
direction. However, the experimental path in the parameter space as a function of 
E is not known and not all of the different scaling regimes correspond to actual 
physical situations. The scaling results found by the RG analysis are verified by 
direct integration of equations, as mentioned in the earlier lectures. A more detailed 
discussion of the analysis and results can be found in our earlier work 40 . 

In constructing equations (60) , we only allowed for local effects, and ignored the 
nonlocalities that are the hallmark of hydrodynamics. One consequence of hydro- 
dynamic interactions is the back- flow velocity in Eq.(58) that can be added to the 
evolution equations (60). Dimensional analysis gives the recursion relation 




(60) 



( m (x,t) m (x',t')) = 2T\\8{x - x')8(t - t'), 
(V±i(x,t)rj ±j (x', t')) = 2T ± 6ij6(x - x')8(t - t'). 



(61) 




<9 7 



7 [zA2-l-(d-2)z/]+0( 7 2 ), 



(62) 
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which implies that, at the nonlinear fixed point, this additional term is surprisingly 
irrelevant for d > 3, and z = 3 due to the nonlinearities. For d < 3, z = d due 
to hydrodynamics, and the nonlinear terms are irrelevant. The situation in three 
dimensions is unclear, but a change in the exponents is unlikely. Similarly, one 
could consider the effect of self-avoidance by including the force generated by a softly 
repulsive contact potential 



The relevance of this term is also controlled by the scaling dimension yf, = vz — \ — [d— 
2)z/, and therefore this effect is marginal in three dimensions at the nonlinear fixed 
point, in contrast with both Rouse and Zimm models where self-avoidance becomes 
relevant below four dimensions. Unfortunately, one is ultimately forced to consider 
non-local and nonlinear terms based on similar grounds, and such terms are indeed 
relevant below four dimensions. In some cases, local or global arclength conservation 
may be an important consideration in writing down a dynamics for the system. 
However, a local description is likely to be more correct in a more complicated system 
with screening effects (motion in a gel that screens hydrodynamic interactions) where 
a first principles approach becomes even more intractable. Therefore, this model is 
an important starting point towards understanding the scaling behavior of polymers 
under a uniform drift, a problem with great technological importance. 
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